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ABSTRACT 

VLBI is a new satellite tracking technique that to date has been applied 
largely to the tracking of synchronous orbits. These orbits are favorable for 
VLBI in that the remote satellite range allows continuous viewing from widely 
separated stations. The primary observable, geometric time delay, is the time 
difference for signal propagation between satellite and baseline terminals. 
Extraordinary accuracy in angular position data on the satellite can be obtained 
by observation from baselines of continental dimensions. In satellite tracking 
though the common objective is to derive orbital elements. A question arises 
as to how the baseline vector bears on the accuracy of determining the elements. 
Our approach to this question is to derive an analytic expression for the time 
delay function in terms of Kepler elements and station coordinates. The 
analysis, which is for simplicity based on elliptic motion, shows that the 
resolution for the inclination of the orbital plane depends on the magnitude of 
the baseline polar component and the resolution for in-plane elements (semi- 
major axis, mean orbital longitude and eccentricity) depends on the magnitude 
of a projected equatorial baseline component. An application is made to obtain 
orbital elements from VLBI tracking of the ATS-3 synchronous satellite. 
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THE VLBI TIME DELAY FUNCTION 
FOR SYNCHRONOUS ORBITS 


INTRODUCTION 

Very Long Baseline Interferometry (VLBI) is a recently developed RF tech- 
nique that is coming into use for the tracking of spacecraft. In this technique 
the spacecraft radiates noise over a broad frequency band and the amplitude 
variations of the signals are simultaneously recorded at two widely separated 
observing stations. Their incoming signals are essentially the same except for 
a difference in arrival time (geometric time delay, r ) and a relative doppler 
frequency shift (fringe rate, ). The stations are each instrumented with VLBI 
backend systems which include independent atomic frequency and timing standards 
and magnetic tape recorders for the processing and recording of the signals. 

A detailed discussion of the methods of VLBI whereby signals are processed and 
reduced to tracking data is given in Ref. 1. Here we are concerned with the 
analytic properties of the VLBI tracking data as related to satellite orbital 
elements. 

The time delay measures the time difference for signal propagation between 
satellite and the terminals of the baseline formed by the two stations. Thus 


Rj and R 2 denoting the range to the stations and c the speed of light. Fringe 
rate is defined as a delay rate. The equation is 

V =ff g (2) 

f being the signal frequency. 

When tracking at remote range these observables essentially yield angular-^ 
data. The resolution for determining position will depend not only on instrumental 
factors such as the stability of the atomic standards, signal recording bandwidth, 
and SNR, but also on the station locations. An electrical interconnection between 
station terminals is not utilized so the baseline can be of continental dimensions. 
It is then possible to achieve extraordinary accuracy in the measurement of 
angles. This is the well known feature of VLBI. 


1 



For satellite tracking there is a further interest to see how the baseline 
vector, both in magnitude and orientation, appear in the problem of determining 
orbital elements. The orbits tracked by VLBI to the present time are largely of the 
near-geostationary type. This report then derives for these orbits an analytical 
expression for r g in terms of station coordinates and Kepler elements and 
examines its dependence on the baseline components. The analysis assumes 
elliptic motion and is limited to arc length such that secular drift can be treated 
as a small quantity. An application of the analysis is made to the tracking data 
of the ATS- 3 satellite . 


ANALYSIS 

Series expansion of t . Let r denote the satellite position relative to the 
Earth center of mass, O, and Sj and S 2 denote the location of the baseline 
terminals. Then by Eq. 1 



The parameter for the series expansion of the radicals is a ratio on the order of 
the Earth radius to the satellite geocentric radial distance. A series with a 
stronger convergence can be formed by utilizing the property of the synchronous 
orbits. This is accomplished by taking as a point of reference a geostationary 
position near the satellite. 


Let t be a time for the satellite crossing of the right ascending node and 
let a denote an Earth-fixed point P in the equatorial plane such that OP equals 
the semi-major axis (Figure 1). The direction of a is taken to coincide with the 
right ascending node at t = t a . Then at this instant a and the satellite position r 
are colinear. 

The separation vector of the satellite from P is 


and that of station locations similarly are 
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Figure 1. Satellite and Station Vector Geometry. 



Figure 2. Geometry of Orbital Elements. 


3 









CT i = S i - a 


a 2 = S 2 - a 


(5b) 


In terms of these vectors the time delay is rewritten as 


:-R 



(6) 


Then using the expansion 


2 -n 3 


N <1 


2 8 16 


( 7 ) 


Eq. 6 becomes 


cr g =a i -P- -^ 2 ) + 1 p 2 ' 1 1 

2 \ a i ^2 


1 ((P-^) 2 _(P-Z 2 ) 2 ^ plfCp-VJ (p-3r 2 y 
v, <y~ 


fl(L-± 
8 W CT 2 


( 8 ) 


The leading members of this series are to be developed in terms of orbital 
elements. 


Notation for the orbital elements is as follows: 

a semi-major axis 

e eccentricity 

i inclination 

co argument of perigee 

fl right ascension of the ascending node 

M mean anomaly 
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f 

E 

n 


true anomaly 
eccentric anomaly 
mean motion. 


Reference frame. In the Earth-centered inertial reference frame (Figure 2) 
the satellite Cartesian coordinates are given by 




l P 1 P 2 P 3 \ 


r cos f 

X 2 

= 

Qi Q 2 Q 3 


1 • f 
r sin f 

X 3j 


i ^2 ^3 j 


0 , 


( 9 ) 


where 

Pj = cos ft cos a) - cos i sin fi sin oj 
P 2 = - cos ft sin eo - cos i sin ft cos co 
P 3 = sin i sin ft 

Q } = sin ft cos o) + cos i cos ft sin ^ 

Q 2 = - sin ft sin a> + cos i cos ft cos u 

( 10 ) 

Q 3 = - sin i cos ft 
Rj = sin i sin w 
R 2 = sin i cos ^ 

Rg = cos i 

In order to refer the satellite position to P a transformation is first made to an 
Earth-fixed Cartesian system u, v, z where the positive direction of the u -axis 
is directed through P. The relation is 



s i n 6 
cos 6 
0 



( 11 ) 
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the rotation angle 6 being of the form 

6 - co t + a ( 12 ) 

e 

where t is measured from epoch, w c is the Earth’s rotation rate, and a is a 
parameter to be determined. 

From Eqs. 9-11 we have 


u = r [cos (Cl - 6) cos (a ; + f) - cos i sin (Cl - 9) sin (co + f)] 
v — r [sin (Cl - 6) cos (co + f) + cos i cos (Cl - 9) sin (w + f)] 
z = r sin i sin (w + f) (13) 


At time t the crossing of the nodal line 

u (t a ) = r. v (t a ) = z (t a ) = 0 
For a positive crossing the conditions are 


(14) 


f (t ) = — ct>, a - Cl - co t 

v a / e a 


(15) 


Then Eq. 12 becomes 


6 = w e (t - t a ) + Cl 


(16) 


The longitude of the u -axis is clearly 

L = n-a g (t a ) ( ] 

a g denoting the right ascension of the Greenwich meridian. 

To determine t we need the relation between the mean anomaly and true 

a 

anomaly which to first order in the eccentricity is 


or 


M = f - 2 e sin f 


(18a) 


n t + M„ = f - 2 e sin f 


(18b) 
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where Mg is the mean anomaly at epoch. Evaluating this at t = t g and using 
Eq. 15 gives the desired relation for t a 

n t a + Mq = - w + 2 e s i n oo 
Now introducing Eq. 19 into 17 yields 

L^n+w + Mg + nt a -a g (t a )-2esinw 

or 

L = fl + w + M 0 -a g (0) - A 

where 

A = (oj -rV)t + 2 e sin w 

' e ' : 8 

In Eq. 20b L is expressed in terms of Q + co+ M, known as the mean orbital 
longitude, and a g both at epoch; and A which is a higher order term. 

The parameter L is a reference longitude for the satellite position. In the 
case of a geostationary orbit it is the one parameter determined from time delay 
measurements. 

Separation vector, p. The near-geostationary orbits which we analyze are 
characterized by the conditions, e< <1, i<<l, and [ ( — n) 1 1 < < | . Restrictions 

on these parameters, form a limit on the magnitude of p„ 


(20a) 

(20b) 

( 21 ) 


By Eq. 5a 



l'«l 


u - a 





p = 

P V 


V 


kl 


1' z 


( 22 ) 


The components of p are to be expressed in terms of the mean anomaly. But 
first the relation between the true and eccentric anomaly 

r cos f = a (cos E - e) 


r sin f = a /l - e 2 sin E (23) 
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is introduced in Eq. 9. Then the components of p in Eq. 22 can be written as 

P u /a= (cos E - e) [cos (fi + co - 9) + (1 - cos i) sin (fi - £) sinoi] 

- /l - sin E [sin (Cl + co - 9) - ( 1 - cos i) sin (Cl - 9) cos co] - 1 


P v /& = (cos E - e) [s in (Cl + co - 9) - (1 - cos i) cos (Cl -9) sin co] 
+ /l - sin E [cos (Cl + co - 9) - (1 - cos i) cos (Cl-9) cos co] 


p x / a = sin i [(cos E - e) sin oj + /l - e 2 sin E cos co] (24) 

The magnitude squared of p appears in one of the terms in the expansion 
of given in Eq. 8. This is readily calculated from Eq. 24 giving 

p 2 / a 2 = 1 + (1 - e cos E) 2 

- 2 {(cos E - e) [cos (0 + co - 9) + (1 - cos i) sin (fi - 6) sin a] 

- /l - sin E [sin (Cl + co - 9) - (1 - cos i) sin (Cl - 9 ) cos a>]} (25) 

The mean anomaly is introduced by the Fourier series expansion of cos E 
and sin E 


cos E - e 


3 e 




cos 2 M 


3 e 2 

+ 

8 


cos 3 M + 


(26a) 


sin E = 1 - 


sin M + ( — - 

2 6 


3 e 2 

sin 2 M + sin 3 M + ■ 


(26b) 


and 
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(26c) 



Combining Eqs. 25, 26a, and 26b and using the approximation, a < < 1 we have 
after rearranging terms 


= 1-cos(M + 0 + oj-(9) 

2 a 2 


+ e 


3 1 

- cos M + — cos (fi + co - 9) - _ cos (2 M + 0 + co - 9) 


13 3 

- - + cos 2 M - — cos (M + D. + co - 6) + —cos (3 M + fi + co - 9) 


— sin M sin (fi + oJ-0) 
2 


— — sin (M + oS) sin (H - 6) + 0 (e 3 , i 2 e) 


(27) 


This equation appears to have terms of the zeroth and first order. A closer 
inspection, however, reveals that the lower order terms cancel. To show this 
we need to rewrite the form of the arguments. The relation used is obtained by 
combining Eqs. 16 and 19 giving 


pM + n + w- <9 = A-(w e -pn)t + (p-l)M 0 (28a) 

or alternatively 

p M + Q. + oj - 6 - (p - 1) M + A - (oj b - n) t (28b) 


where p is a positive or negative integer. 
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For the first two terms of Eq. 27 we have using Eq. 28a 

l-cos(M + fi + oj-(9)- i - cos [A - (co e - n) t] 

. [A - (co - n) t] 2 


(29) 


This is a second order term since we consider arc length such that | (w e - n) 1 1 <<1. 
Next the coefficient of e in Eq. 27 is 


^ 1 

I = - cos M + — cos (Q. + co - 6) - —cos + + 

2 2 


(30a) 


or rearranging the arguments by Eq. 28b 


3 1 

I=-cosM + — cos [M+(w e -n)t-A] cos [M - (w e - n). t + A] 

2 2 


Then using the approximation cos ( M + 77 ) = cos M - 77 sin M we have 


I = - 2 [(co - n) (t - t ) - 2 e sin a j] sin M 


(30b) 


Thus by a cancellation of the zeroth order terms the expression of Eq. 30a 
becomes of the first order. 

We now rewrite Eq. 27 using Eq. 28b 


p 2 

= 1 - cos [A - (co - n) t] + e 


2 a 2 


-j- co 


s M + — COS [M + (ct> e - n) t - A] 


- —cos [M - (« g - n) t + All- + — — i- — - — cos [A - (co - n) t] 
* J 2 l 2 4 


3 1 

+ cos 2 M + — cos [2 M + A - (w e - n) t] - — cos [A - (co g - n) t] 
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+ —cos [2 M - A + (C0 e - n) t] 


- cos [2 M + 2 a) 


H{ 

-A + (a) e -n)t]| + 0 


cos [A - (co e - n) t] 
(e 3 , i 2 e) 


( 31 ) 


bearing in mind that the coefficient of e can be approximated by the expression 
of Eq. 30b. 

Thus p 2 is second order in e, i, and (a) e - n)t. Now the term in the expansion 
of r g (Eq. 8) 


a l ^2 


is still another order of magnitude lower since a l = a 2 . It is generally small 
compared to the linear term, p • (a l - <x 2 ). To analyze the latter we need to 
consider the vector difference, & = 5- 1 - 5- 


Sigma vectors. By Eq. 5b 


M 


1 S x jCosL + S yl sinL-a 

^.v 

= 

- S x i sin L + S y j cosL 



1 s --‘ i 


and 


(32a) 


1 

j a 2.u| 


1 S x, 2 cos L + S y >2 sin L - a 

ii 

°2, v 

= 

- S x, 2 Sin L + S y,2 COS L 


[ Cr I 


1 S „ 


i 2 ’*l 


l z - 2 


(32b) 


l * ®y. l » ® *, i an< ^ ®x, 2 > ® y , 2 * 2 denoting the Earth-centered Cartesian 

station coordinates. The magnitudes are 


°i = /a 2 + Sjf - 2 a (S x>1 


cos L + S j sin L) 


(33a) 
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= /a2 + s j - 2 a CS x 2 cos L + S y 2 sin L) 


(33b) 


where 


s . = »SJ^T^, I + sJ. 1 




The difference of the unit vectors is denoted by 


a i ~ a 2 = = 


(34) 


the components being 


cr u = S x cos L + S y sinL - a 
cr v = -S sinL + S cosL 


(35) 


where 


and 


S 


X , 1 X , 2 


cr. 


c _ S y-i S y.2 ^ s *,i S z, : 

’ i O, 


CT 1 ^2 


a = a 


°1 C7 2 


(36) 


The magnitude of a is then 


<r = /a 2 + S* + S* + S* - 2 a (S x cos L + S sinL) 


(37) 
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It will be useful to express the components of a through angles /i, v defined 
by the equations 


o" u = cr cos yLt cos v 

cr v = ct c os /i sin y (38) 


= cr sin /i 


Then 


CT 

t an v - 

a 

u 

and 


(39) 


ct 

■sin fi = __L 


CT 


(40) 


The leading term in the expansion of in Eq. 8 is ct^ - o- its differential 

is readily determined by the above equations to be 


A (CX l ~ a 2 ) = “ a u A 3 ~ 3 °v A L ( 41 ) 

The geometrical significance of the components of c? is clarified when 

°i = a 2 * This is a reasona ble approximation when the angle between a and 
the station directions is about 50° or less. Then 



But by Eq. 5 b 



where D is the baseline vector. Hence 


cr cr = D 

1 u u 


o- cr -= D 

1 V V 


a i a z = D. 


(42) 


(5c) 


(43) 
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and 


• D v . D z D 

tan v = - , sin /i = — , cos u = — — L 

D D D 


(44) 


where D eq = / + D ^ is the projection of D on the equatorial plane. 


Finally another useful relation is 


a cos fi = 


D COS fjL _ 


(45) 


Expansion of p • o- . The term of r (Eq. 8) linear in p is of the form 
P' (3- - a )= p • a. From Eqs, 24 and 38 we have that 

p - a - a a cos p. {- cos v + (cos E - e) (cos (fi + co - 9 - v) 


+ (1 - cos i) sin (CL - 6 - v) sinw) - /l - e 2 sinE(sin(fi + oj-i9-v) 


- (1 - cos i) sin (fi - 9 - v) cos u)} 


+ a cr sin/i sin i [(cos E - e) sin w + /l - si n E cos oo) 


(46) 


Then eliminating E by Eqs. 24a and 24c and using the approximation i < < 1 


:os p. 


p • a - aacos/i(- cosv + cos (F + M) ■‘r ~ (cos (F - M - 2 co) 


3 1 

- cos (F + M)) + e (-— cos F + j cos (F + 2 M) 


e i 


+ — [cos (F - 2 M - 2 co) - cos (F + 2 M) - 3 cos (F - 2 co) 


+ 3 cos F)] + — - [- 4 cos (F + M) + 3 cos (F + 3 M) + cos (F - M)] 


') 


(47) 
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a a s 1 


f 3 e / 3 e 2 \ 

inpsiniV- sin w + [ 1 - 1 sin 

l 2 \ 8 / 


(M + co) 


e . e 2 3 e 2 

+ — sin (2 M + co) sin M cos co + sin (3M + co) 

2 4 8 


} 


where 


F - Q. + co-9-v (48) 

To show the dependence of the arguments in Eq. 47 on the mean anomaly, Eq. 48 
can be rewritten using Eq. 28b 

F + p M = (p - 1) M + A - (co e - n) t ( 49 > 

where the right member is now a sum of an integral multiple of M plus higher 
order terms. 

It is useful to write p •<? in a form that shows the harmonic of each term 
explicitly. From Eq. 49 we readily have 

F + p M = - O e - p n) t + A + (p - 1) ( 50 ) 

Introducing this into Eq. 47, yields 

p:a - a a cos p cos v. + cos (A - v) cos (co g - n) t + sin (A - v) sin (w e -n) t 
+ y j- 3cos(A-^-)^)cosw e t - 3 sin (A - y - M Q ) sin co e t 


+ cos (A - v + M Q ) cos (w e - 2 n) t + sin (A - v + M Q ) sin ( co g 




- 4 cos (A - i^) cos (<x> e 


- n) t - 4 sin (A - v) sin (a> e 


- n) t 


+ 3 cos (A - v + 2 M 0 ) cos (co e - 3 n) t + 3 sin (A - v + 2 M 0 ) sin ( u> e - 3n)t 
+ cos (A - V - 2 M 0 ) cos (oo e + n) t + s in (A - v - 2 M ff ) sin (co e + n) t"J 
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i 2 r 

+ “ l cos (A - v - 2 M„ - 2 co) cos (o> e + n) t + sin (A - v - 2 M Q - 2 co) sin (co g + n) t 
- cos (A - V) cos (co e - n) t - sin (A - v) sin (co e - n) tj 


e 2 i 
8 


3 cos (A - v - M 0 ) cosw e t + 3 sin (A - v - M Q ) sin w e t 


- 3 cos (A - v - M 0 - 2 a) cos t - s in (A - v - M q - 2 oS) sin a> e t 

+ cos (A - v - 3 M 0 - 2 oS) cos (o) e + 2 n) t + sin (A - v - 3 M Q -2 oS) sin (c<J e+ 2n) t 

- cos (A - v + M Q ) cos (co e - 2 n) t - sin (A - v + M Q ) sin (o> e - 2 n) t 
+ a a sin n sin i jsin (co + M 0 ) cos n t + cos (co + M Q ) sin n t 


£ 

2 


+ - - 3 sin co + sin (co + 2 M Q ) cos 2 n t + cos (co + 2 M Q ) sin 2 n t 


[3 . 

-si 

L 2 


— —sin (a. + M 0 ) cos n t + - cos (co + M 0 ) sin n t 


3 . 3 

_ 2 - sln + 3 M„) cos 3 n t - - cos (co + 3 M 0 ) sin 3 n t 


+ cos a) cos M 0 sin n t + cos co sin M Q cos 


nt ]} 


(51) 


Discussion . The above equation is written as the sum of two parts, one 
having the common factor a o- cos y. and the other having the factor a a- sin fi sin i . 
Within the approximation 0- = a 2 it can, by Eqs. 44, 45, and 51, be put in 
the form 


p-v = — D, 


cr. 


eq 



+ — D sin i 


(52) 
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the ratio a/ o r 1 ~ 0.9. Thus the members on the right side depend on the equatorial 
component of the baseline, D eq , and the polar component, D z , respectively. This 
expression affords some incite on how these components bear on determining 
orbital elements. 

Consider first r g measured from a polar oriented baseline. With reference 
to Eqs. 51 and 52 we can see that the resolution for determining the inclination 
depends on the baseline according to the formula a i ~ (D z cos i)" 1 . This is an 
analogue for the orbital element of the more familiar relation in VLBI where 
resolution in measuring angular position depends on the reciprocal projected 
baseline length. The determination of elements other than the inclination varies 
as (D z sin i) -1 which is much degraded for small inclination. The 
polar oriented baseline is then primarily for accurate measurement of the inclina- 
tion of the orbital plane. 


When measured from an equatorial oriented baseline r has a term linear 
in time. By Eq. 44 this term can be rewritten 


a a cos /x sin (A - v) sin (aj e - n) t 


= ~ D eq [sin (A - v)] ( co e - n) t 
= — D (co - n) t (53) 

Or v e 

since- | A | < < 1 and | (« e - n)t | < < 1. The difference o> e - n then can be measured 
as D” 1 , that is the reciprocal of the projected equatorial baseline component normal 
to the OP line. This term is the basis for a precise determination of n since it 
appears differenced with « e . Further examination of the terms of Eq. 51 cor- 
responding to D eq shows that the eccentricity is of the first order and inclination 
is of the second order. 

The elements Cl , u > , and M 0 arise in the arguments of trigonometric terms 
having coefficients e, i, or (w, - n) t to the first order. , Consequently their 
resolution is degraded. For the limiting case of a geostationary orbit these 
elements become indeterminate. 

The mean orbital longitude is however well defined. This appears in the 
leading term of the expansion of r g , namely a - 1 - o 2 . By Eqs. 41 and 43 
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A a - 


(54) 


A (o-j - a 2 ) = - ■ 


D 


a D 


•A (fi + co + M„) 


The semi-major axis can be considered determined from the mean motion as 
described above. Then 


A (fi + & + M 0 ) ~ - A (o-j -a 2 )/D v 

It is possible then to determine the summed elements, fi + co + M 0 , with a higher 
precision than its individual components. 


TRACKING OF ATS-3 

VLBI satellite tracking observations yield geometric time delay data as a 
function of time. By identifying this numerical function with the analytical 
expression, t (t) , the Kepler orbital elements can be deduced. This has been 
done using data from a VLBI tracking experiment on the ATS-3 communication 
satellite. 

The C-band noise signals of the spacecraft were observed at the ATS sites 
of Rosman, North Carolina and Mojave, California during June 3, 1971 to 
June 10, 1971 when the satellite longitude was approximately 79°W, a few 
degrees east of Rosman. The difference in the latitude of the two stations 
differs by scarcely more than 8', the geodetic coordinates being 


Latitude 


Longitude (W) 


Height (m) 


Mojave 35 0 19’53V78 116°53 , 16'.'14 

Rosman 35°11 , 56'.'07 82 0 52'32’.'81 


887 

828 


The equatorial component is 3050 km and polar component is 12 km. 

Figure 3 shows the distribution of tracking periods superimposed on a 
curve representing the time delay function. The periods are each of three 
minute duration. On day numbers 152, 153, 154, and 155 observations were for 
a restricted part of the day but for day 159 observations were distributed over 
a twelve hour period. The precision of the time delay data is 5 to 10 nanoseconds 
corresponding to 1.5 to 3 meters. 
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•TRACKING PERIOD 



Figure 3. 


Time delay curve showing the VLBI tracking periods of ATS-3, days 


152-159, 1971. 


As a starting point for calculating the orbital elements the terms of the time 
delay function were evaluated using a priori elements (based on range and range 
rate tracking from Mojave and Rosman) issued by the ATS project. Terms 
greater than 30 nanoseconds (10 meters) were retained for which the functional 
form was essentially 

cr = c + c t+c„ cos co t + c. s in co t + c cos 2 co t 

g 0 1 ^ c 

+ c s in 2 co t (55) 

5 e 

The coefficients were evaluated by the method of least squares, the standard 
deviation in time delay being 0.66 microseconds. One of the coefficients was 
smaller than the st. dev. and therefore was not regarded as significant. Five 
elements were determined and are compared with those issued by the ATS project. 
The epoch is May 31, 1971, 00 h . 


Element 

ATS Project 

VLBI Time Delay 

Difference 

a (km) 

42,165.43 

42,165.54 

0.11 

e 

2.914 X 10 -3 

2.96 x Id' 3 

5 x 10' 5 

i [ 

If 706 

If 6 

O 
• o 
H- 1 

M o 

87° 150 

95° 

8°' 

1 + CO + M„ 

16 8 f 65 2 

168 °61 

-Of 04 

0 

CO 

-0?712 

— 

— 
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It should be understood that the method here is not expected to yield ac- 
curate values for the elements since the model is based on the approximation of 
elliptic motion. To derive a best estimate of the ATS-3 orbit from VLBI obser- 
vations we are using GEODYN (Ref. 1). This program has an optimum force 
field model and employs rigorous numerical methods which are designed for 
application to geodynamic investigations. 


ACKNOWLEDGMENTS 

The author gratefully acknowledges many interesting discussions on this 
work with Drs. J. Ramasastry and S. Ross. 


REFERENCE 

1. Ramasastry, J., et al., Tracking of the ATS-3 Synchronous Satellite by 
the Very Long Baseline Interferometer (VLBI) Technique, GSFC 
X-553-72-290, July 1972. 


nasa-gsfc 


20 



